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Abstract. We derive the value of Hq using the inverse diame- 
ter and magnitude B-band TuUy-Fisher relations and the large 
all-sky sample KLUN (5171 spiral galaxies). Our kinematical 
model was that of Peebles centered at Virgo. Our calibrator 
sample consisted of 15 field galaxies with cepheid distance 
moduli measured mostly with HST. A straightforward appli- 
cation of the inverse relation yielded Hq w 80 kms^^ Mpc^^ 
for the diameter relation and Hq ~ 70 kms^^ Mpc~^ for the 
magnitude relation. Hq from diameters is about 50 percent and 
from magnitudes about 30 percent larger than the correspond- 
ing direct estimates (cf. Theureau et al. 1997b). This discrep- 



ancy could not be resolved in terms of a selection effect in 
log Knax nor by the dependence of the zero-point on the Hub- 
ble type. 

We sh owed that a new, calibrator selection bias (Teeriko- 
rpi 1999), is present. By using samples of signicificant size 
(N=2142 for diameters and N=1713 for magnitudes) we found 
for a homogeneous distribution of galaxies (a — 0): 



B- 



- Hq = 52_4 km Mpc^ for the inverse diameter 
band TuUy-Fisher relation, and 

- Hq = 53^5 kms^^ Mpc^^ for the inverse magnitude B 
band TuUy-Fisher relation. 



Also Hq's from a fractal distribution of galaxies (decreasing 
radial number density gradient a — 0.8) agree with the direct 
predictions. This is the first time when the inverse TuUy-Fisher 
relation clearly lends credence to small values of the Hubble 
constant Hq and to long cosmological distance scale consis- 
tently supported by Sandage and his collaborators. 

Key words: Galaxies: spiral - Galaxies: distances and redshifts 
- Cosmology: distance scale 



1. Introduction 

The determination of the value of the Hubble constant, Hq, is 
one of the classical tasks of observational cosmology. In the 
framework of the expanding space paradigm it provides a mea- 
sure of the distance scale in FRW universes and its recipro- 
cal gives the time scale. This problem has been approached 
in various ways. A review on the recent determinations of 
the value of Hq shows that most methods provide values at 
Hq ~ 55 ... 75 (for brevity we omit the units; all Hq val- 
ues are in kms^^ Mpc^^): Virgo cluster yields 55 ± 7 and 
clusters from Hubble diagram with relative distances to Virgo 
57 ± 7 (Federspiel et al. 1998 ), type la supernovae give 60 ± 10 
(Branch |1998| ) or 65 ± 7 (Riess et al. |1998| ), Tully- Fisher rela- 
tion in I-band yields 69 ± 5 (Giovanelli et al. 1997 ) and 55 ± 7 
in B-band (Theureau et al. 



1997b, value and errors combined 



from the diameter and magnitude relations), red giant branch 
tip gives 60 ± 11 (Salaris & Cassisi 1998), gravitational lens 
time delays 64± 13 (Kundic et al. 1997) and the 'sosies' galaxy 
method 60 ± 10 (Paturel et al. |l998t ). Sunyaev -Zeld ovich effect 
has given lower values , 49 ± 29 by Cooray ( |1998[ ), 47^11 by 
Hughes & Birkinshaw ( 1998 1, but the uncertainties in these re- 



sults are large due to various systematical effects (Cen, 1998) 



Surface brightness fluctuation studies provide a higher value of 



87 ± 11 (Jensen et al., 1999 1, but most methods seem to fit in 
the range 55 - 75 stated above. An important comparison to 
these local values may be found after the cosmic microwave 
background anisotropy probes (MAP and Planck) and galaxy 
redshift surveys (2dF and SPSS) offer us a multitude of high 
resolution data (Eisenstein et al., 1998). Note that most of the 
errors cited here as well as given in the present paper are la 
errors. 

The present line of research has its roots in the work of 
BottinelU et al. (1986), where Hq was determined using spiral 
galaxies in the fie Id. They used the direct Tully-Fisher relation 
(Tully & Fisher [1977|): 
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where M is the absolute magnitude in a given band and 
log Vinax is the maximum rotational velocity measured from 
the hydrogen 21 cm line width of each galaxy. Gouguenheim 



( 1969) was the first to suggest that such a relation might exist 



as a distance indicator. 



ination of the so-called Malmquist bias. In general terms, the 
determination of Hq is subject to the Malmquist bias of the 2"*^ 
kind: the inferred value of Hq depends on the distributio n of th e 
derived distances r for each true distance r' (Teerikorpi 1997 ). 
Consider the expectation value of the derived distance r at a 
given true distance r': 



E{r\r') = j drrP{r\r'). 



(2) 



The integral is done over derived distances r. For example, con- 
sider a strict magnitude limit: for each true distance the derived 
distances are exposed to an upper cut-off. Hence the expecta- 
tion value for the derived distance r at r' is too small and thus 
Hq will be overestimated. 

Observationally, the direct Tully-Fisher relation takes the 
form: 



X = slope X p + est, 



(3) 



where we have adopted a shorthand p for log Knax and X de- 
notes either the absolute magnitude M or log D, where D la- 
bels the absolute linear size of a galaxy in kpc. In the direct 
approach the slope is determined from the linear regression of 
X against p. The resulting direct Tully-Fisher relation can be 
expressed as 



E{X\p) = ap + b. 



(4) 



Consider now the observed average of X at each p, {X)p, as a 
function of the true distance. The limit in x (the observational 
counterpart of X) cuts off progressively more and more of the 
distribution function of X for a constant p. Assuming X = 
log D one finds: 



{X)p > EiX\p), 



(5) 



The inequality gives a practical measure of the Malmquist bias 
depending primarily on p, r', ax and xu^. The equality holds 
only when the x-limit cuts the luminosity function $(X) in- 
significantly. 

That the direct relation is inevitably biased by its nature 
forces one either to look for an unbiased subsample or to find 
an appropriate correction for the bias. The former was the strat- 
egy chosen by Bottinelli et al. (1986) where the method of 



normalized distances was introduced. This is the method cho- 
sen also by the KLUN project. KLUN (Kinematics of the Loal 
Universe) is based on a large sample, which consists of 5171 
galaxies of Hubble types T=l-8 distributed on the whole celes- 
tial sphere (cf e.g. P a turel|1994 Theureau et al. |1997b[ ). 



Sandage (1994a 1994b) has also studied the latter ap- 



only on the imposed x-limit but also on the rotational velocities 
and distances, he introduced the triple-entry correction method, 
which has consistently predicted values of Hq supporting the 
long cosmological distance scale. As a practical example of 
this approach to the Malmquist bias cf. e.g. Federspiel et al. 



Bottinelli et al. ( 1986 ) paid particular attention to the elim- ( 1994 ) 



Bottinein et al. (1986) found Hq = 72kms"^ Mpe"^ us- 



ing the method of normalized distances, i.e. using a sample 
cleaned of galaxies suffering from the Malmquist bias. This 
value was based on the de Vaucouleurs calibrator distances. 
If, instead, the Sandage-Tammann calibrator distances were 



used BottinelH et al. (1986) found Hq = eSkms^^Mpe- 



(or Hq ~ 56kms^^ Mpc^^ if using the old ST calibration). 
One appreciates the debilitating effect of the Malmquist bias 
by noting that when it is ignored the de Vaucouleurs calibra- 
tion yields much lar ger val ues: Hq ~ lOOkms"^ Mpe^^. 

Theureau et al. (1997b) by following the guidelines set out 
by Bottinelli et al. (1986) determined the value of Hq using the 
KLUN sample. Hq was determined not only using magnitudes 
but also diameters because the KLUN sample is constructed to 
be complete in angular diameters rather than magnitudes (com- 
pleteness limit is estimated to be D25 — 1.6). Left with 400 
unbiased galaxies (about ten times more than Bottinelli et al. 
(|198^ were able to use) reaching up to 2000-3000 km s~^ 



they found using the most recent calibration based on HST ob- 
servations of extragalactic cepheids 

- Hq = 53.4 ± 5.0 km s^^ Mpe^^ from the magnitude rela- 
tion, and 

- Hq = 56.7±4.9kms^^ Mpe^^ fromthediameterrelation. 

They also discussed in their Sect. 4.2 how these results change 
if the older calibrations were used. For example, the de Vau- 
couleurs calibration would increase these values by 11 %. We 
expect that a similar effect would be observed also in the 
present case. 

In the present paper we ask whether the results of Theureau 
et al. (1997b) could be confirmed by implementing the inverse 
Tully-Fisher relation: 



p = a X + h\ 



(6) 



This problem has special importance because of the "unbiased" 
nature that has often been ascribed to the inverse Tully-Fisher 
relation as a distance indicator and because of the large number 
of galaxies available contrary to the direct approach where one 
is constrained to the so call ed unbi ased plateau (cf. Bottinelli 
et al. 1986; Theureau et al. 1997b). The fact that the inverse 



relation has it own particular biases has recei ved increasi ng at 
tention during the years ( Fouqu e et al 
Willick 



1991 



19901 , Teerikorpi 
Teerikorpi 1993| , Ekholm & Teerikorpi 



Freudling et al. 1995, Ekholm & Teerikorpi 1997, Teerikorpi 



1990 



1994 



et al. 1999 and, of course, the present paper). 



proach. By recognizing that the Malmquist bias depends not 



2. Outlining the approach 

As noted in the introduction the KLUN project approaches the 
problem of the determination of the value of Hq using field 
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galaxies with photometric distances. Such an approach reduces 
to three steps 

1 . construction of a relative kinematical distance scale, 

2. construction of a relative redshift-independent distance 
scale, and 

3. establishment of an absolute calibration. 

Below we comment on the first two steps. In particular we fur- 
ther develop the concept of a relevant inverse slope which may 
differ from the theoretical slope, but is still the slope to be used. 
The third step is addressed in Sect. 6. It is hoped that this re- 
view clarifies the methodological basis of the KLUN project 
and also makes the notation used more familiar. 

2.1. The kinematical distance scale 

The first step takes its simplest form by assuming the strictly 
linear Hubble law; 



^kin — Vo/Hq 



(7) 



where Vo is the radial velocity inferred from the observed red- 
shifts and Hq is some input value for the Hubble constant. Be- 
cause Vo reflects the true kinematical distance R^^^^ via the true 
Hubble constant HX 



Vo/H*, (8) 
one recognizes that Eq. ^ sets up a relative distance scale: 



Rkii 

R 



'kin 



HI 



(9) 



In other words, log dkin is known next to a constant. 

In a more realistic case one ought to consider also the pecu- 
liar velocity field. In KLUN one assumes that pecuhar veloci- 
ties are governed mainly by the Virgo supercluster. 

In KLUN the kinematical distances are inferred from Vq's 
by implementing the spherically symmetric model of Peebles 
(1976) valid in the linear regime. In the adopted form of this 
model (for the equations to be solved cf. e.g. Bottinelli et al. 



1986, Ekholm 1996) the centre of the peculiar velocity field 



is marked by the pair of giant ellipticals M86/87 positioned at 
some unknown true distance R* which is used to normalize the 
kinematical distance scale: the centre is at a distance dkin = 1- 
The required cosmological velocities V^oi (observed veloc- 
ities corrected for peculiar motions) are calculated as 



Vcoi — Ci 



(10) 



where the constant Ci defines the linear recession velocity of 
the centre of the system assumed to be at rest with respect to 
the quiescent Hubble flow: 



(11) 



Vo (Vir) is the presumed velocity of the centre and V^^p is the 
presumed infall velocity of the Local Group into the centre of 
the system. 



2.2. The redshift-independent distances 

The direct TuUy-Fisher relation is quite sensitivite to the sam- 
pling of the luminosity function. On the other hand, when im- 
plementing the inverse Tully-Fisher relation (Eq. ||) under ideal 
conditions it does not matter how we sample X (Schechter 



1980) in order to obtain an unbiased estimate for the inverse pa- 



rameters and, furthermore, the expectation value E{r\r') is also 
unbiased (Teerikorpi 1984). However, we should sample all 
log Knax/c" each constant true X in the sample. This theoreti- 
cal prerequisition is often tacitly assumed in practice. For more 
formal treatments on the inverse relation cf. Teerikorpi(1984, 



1990, 1997) and e.g. Hendry & Simmons (1994) or Rauzy & 



Triay (|199§). 

In the inverse approach the distance indicator is 

X ^ A'(p)x +cst., 



(12) 



where A' — 1 /a' fo llowing the notation adopted by Ekholm 
& Teerikorpi ( 1997 ; hereafter ET97). The inverse regression 
slope a' is expected to fulfill 



{p)x = E{p\X) ^ a' X + est. 



(13) 



{p)x is the observed average p for a given X. Eq. |l3] tells that 
in order to find the correct a' one must sample the distribu- 
tion function 4>x{p) in such a way that {p)x — iPo)x, where 
{po)x is the central value of the underlying distribution func- 
tion. 4>x{j>) is presumed to be symmetric about (po)x for all 
X. ET97 demonstrated how under these ideal conditions the 
derived log Hq as a function of the kinematical distance should 
run horizontally as the adopted slope approaches the ideal, the- 
oretical slope. 

In practice the parameters involved are subject to uncertain- 
ties, in which case one should use instead of the unknown the- 
oretical slope a slope which we call the relevant inverse slope. 
We would like to clarify in accurate terms the meaning of this 
slope which differs from the theoretical slope and which has 
been more heuristically discussed by Teerikorpi et al. (1999). 



The difference between the theoretical and the relevant slope 
can be expressed in the following formal way. Define the ob- 
served parameters as 



Xo = X 



Po=P 



fkir 



(14) 
(15) 



where X is inferred from x with a measurement error ex and 
the kinematical distance dkin has an error ekin due to uncertain- 
ties in the kinematical distance scale, is the observational 
error on p. The theoretical slope a[ is|^ 

, _ Cov{X,p) 



Gov{X, X) ' 



(17) 



' We make use of the formal definition of the slope of the linear 
regression of y against x with 



Cov{x,y) = 



J2{x - {x))iy - jy)) 
(N-l) 



(16) 
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while the observed slope is 

Cov{Xo,Po) Cov(X,p) + Cov(e^ + ekin,ep) 



cr£ + a 



(18) 



kin 



° Cov{Xo,Xo) Cov{X,X) 
We call the slope a'^ relevant if it verifies for all Xo (Eq. |13|) 

{Po)x^ = E{p\X,) = <Xo + est. (19) 



This definition means that the average observed value of po at 
each fixed value of Xo (derived from observations and the kine- 
matical distance scale) is correctly predicted by Eq. IS 
also that in the case of diameter relation, e^, Ckin and e 



Note 
are 



only weakly correlated. Thus the difference between the rele- 
vant slope and the theoretical slope is dominated by cr^ + ^Tkin- 
In the special case where the galaxies are in one cluster (i.e. 
at the same true distance), the dispersion dkin vanishes. In or- 
der to make the relevant slope more tangible we demonstrate in 
Appendix A how it indeed is the one to be used for the deter- 
mination of Hq. 

Finally, also selection in p and type effect may affect the de - 
rived slope making it even shallower. Theureau et al. (1997c) 
showed that a type effect exists seen as degenerate values of 
p for each constant linear diameter X. Early Hubble types ro- 
tate faster than late types. In addition, based on an observa- 
tional program of 2700 galaxies with the Nangay radiotele- 
scope, Theureau (1998) warned that the detection rate in HI 



varies continuosly from early to late types and that on average 
~ 10% of the objects remain unsuccessfully observed. Influ- 
ence of such a selection, which concerns principally the ex- 
treme values of the distribution function (j){p), was discussed 
analytically by Teerikorpi et al. (1999). 



3. A straightforward derivation of log Ho 

3.1. The sample 



KLUN sample is - according to Theureau et al. (1997b) - com- 
plete up to — 13'."25, where is the corrected total B- 
band magnitude and down to logl^js = 1-2, where Z?25 is 
the corrected angular B-band diameter The KLUN sample was 
subjected to exclusion of low-latitude {\b\ > 15°) and face-on 
(log i?25 > 0.07) galaxies. The centre of the spherically sym- 
metric peculiar velocity field was positioned at / — 284° and 
b = 74°. The constant Ci needed in Eq. |lo| for cosmologi- 
cal velocities was chosen to be 1200 km s^^ with 14(Vir) — 
980kms-i and y^G ^ 220 km s^^ (cf. Eq. |ll]). After the 
exclusion of triple-valued solutions to the Peebles' model and 
when the photometric completeness limits cited were imposed 
on the remaining sample one was left with 1713 galaxies for 
the magnitude sample and with 2822 galaxies for the diameter 
sample. 

3.2. The inverse slopes and calibration of zero-points 



Theureau et al. (1997ai derived a common inverse diameter 
slope a' ~ 0.50 and inverse magnitude slope a' k, —0.10 for 
all Hubble types considered i.e. T=l-8. These slopes were also 
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Fig. 1. The slope a' = 0.50 forced to the calibrator sample 
with Cepheid distances yielding b'^^^ = 1.450, when no type 
corrections were made. 
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Fig. 2. The slope a' = —0.10 forced to the calibrator sample 
with Cepheid distances yielding b'^^^ = 0.117, when no type 
corrections were made. 

show n to obe y a a simple mass-luminosity model (cf. Theureau 
et al. 1997a I. With these estimates for the inverse slope the 



relation can be calibrated. At this point of derivation we ig- 
nore the effects of type-dependence and possible selection in 
log Vniax- The calibration was done by forcing the slope to the 
calibrator sample of 15 field galaxies with cepheid distances, 
mostly from the HST programs (Theureau et al. 
Table 1.). The absolute zero-point is given by 



1997b 



cf. their 



Km- 



E(logKnax-a'^) 



(20) 



where the adopted inverse slope a' — 0.50 yields b'^^^ — 1.450 
and a' — —0.10 b'^^^ — 0.117. In Fig. [l| we show the calibra- 
tion for the diameter relation and in Fig. || for the magnitude 
relation. 
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Fig. 3. Panel (a): The log Hq vs. Vcor diagram for the calibrated 
inverse Tully-Fisher relation logVmax = 0.50 log 13 + 1.450. 
The horizontal solid line corresponds to the average value 
(logffo) — 1.92. Panel (b): the average values (logi?o) (cir- 
cles) are shown as well as the average of the whole sample. The 
averages were calculated for velocity bins of size 1000 km s~^. 
Total number of points used was N — 2822. 



3.3. Ho without type corrections 

ET97 discussed in some detail problems which hamper the de- 
termination of the Hubble constant Hq when one applies the in- 
verse Tully-Fisher relation. They concluded that once the rele- 
vant inverse slope is found, the average (log Hq) shows no ten- 
dencies as a function of the distance. Or, in terms of the method 



of normalized distances of Bottinelli et. al. (1986), the unbi- 
ased plateau extends to all distances. ET97 also noted how one 
might simultaneously fine-tune the inverse slope and get an un- 
biased estimate for logiJo. The resulting \ogHo vs. kinemati- 
cal distance diagrams for the inverse diameter relation is given 
in Fig. ^and for the magnitude relation in Fig. ^ Application of 
the parameters given in the previous section yield (logiJo) = 
1.92 correpondingto Hq = 83.2 kms^^ Mpc^^ forthediame- 
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Fig. 4. The sample imposed to the strict magnitude limit = 
13"'25 (N=1713). The forced solution yields (logiJo) = 1-857 
or Ho = 71.9kms-iMpc-i. 



ter sample and (logi?o) = 1-857 or iJo = 71.9kms"i Mpc"^ 
for the magnitude sample. These averages are shown as hori- 
zontal, solid straight lines. In panels (a) individual points are 
plotted and in panels (b) the averages for bins of 1000 kms~^ 
are given as circles. 

Consider first the diameter relation. One clearly sees how 
the average follows a horizontal line up to 9000 km s^^. At 
larger distances, the observed behaviour of (Hq) probably re- 
flects some selection in log Vmax in the sense that there is an 
upper cut-off value for log Vmax- Note also the mild downward 
tendency between 1000 km and 5000 km s^^. Comparison 
of Fig. ^ with Fig. || shows how (log Hq) from magnitudes and 
diameters follow each other quite well as expected (ignoring, of 
course, the vertical shift in the averages). Note how the grow- 
ing tendency of (log Hq) beyond 9000 km s^^ is absent in the 
magnitude sample because of the limiting magnitude: the sam- 
ple is less deep. This suggests that the possible selection bias 
in log V^max does not affect the magnitude sample. 
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One might, by the face-value, be content with the slopes 
adopted as well as with the derived value of Ho. The observed 
behaviour is what ET97 argued to be the prerequisite for an un- 
biased estimate for the Hubble constant: non-horizontal trends 
disappear It is - however - rather disturbing to note that the 
values of Hq obtained via this straightforward application of 
the inverse relation are significantly larger than those reported 



diets some 50 percent larger value and the magnitude relation 
some 30 percent larger value than the corresponding direct re- 
lations. In what follows, we try to understand this discrepancy. 



> 

by Theureau et al. ( 1997b i. The inverse diameter relation pre- o 




4. Is there selection in log 1^ 



0.5 



1 .5 



2.5 



The first explanation coming to mind is that the apparently 
wellbehaving slope a' — 0.5 (a' = —0.1) is incorrect because 
of some selection effect and is thus not relevant in the sense 
discussed in Sect. 2.2 and in Appendix A. The relevant slope 
brings about an unbiased estimate for the Hubble parameter (or 
the Hubble constant if one possesses an ideal calibrator sample) 
if the distribution function of logKnax, 4>{p)x, is completely 
and correctly sampled for each X. Fig. H showed some prelim- 
inary indications that this may not be the case as regards the 
diameter sample. 



Teerikorpi (1999) discussed the effect and significance of a 
strict upper and/or lower cut-off on ip{p)x- For example, an up- 
per cut-off in 4){p)x should yield a too large value of Hq and, 
furthermore, a too shallow slope. Their analytical calculations 
given the gaussianity of 4>{p)x show that this kind of selection 
effect has only a minuscule affect unless the cut-offs are con- 
siderable. Because the selection does not seem to be significant, 
we do not expect much improvement in Hq. 

There is, however, another effect which may alter the slope. 
As mentioned in Sect. 2.2 the type-dependence of the zero- 
point should be taken into account. Because the selection func- 
tion may depend on the morphological type it also affects the 
type corrections. This is clearly seen when one considers how 
the type co rrections are actually calculated. As in Theureau et 
al. (1997b) galaxies are shifted to a common Hubble type 6 by 
applying a correction term Ab' = b'{T) — b'{6) to individual 
log Vmax values, where 



b'{T) = (l0gKnax)T-a'W2 



(21) 



Different morphological types do not have identical spatial oc- 
cupation, which is shown in Fig. ^ for Hubble types 1 and 8 as 
dashed lines corresponding to forced solutions using the com- 
mon slope a' — 0.5. The strict upper and lower cut-offs would 
influence the extreme types more. Hence we must first more 
carefully see if the samples suffer from selection in log VJnax 

The inverse Tully-Fisher diagram for the diameter sample 
is given in Fig. ^ The least squares fit (a' — 0.576, b' = 1.259) 
is shown as a solid line. One finds evidence for both an up- 
per and lower cut-off in the log Vmax-distribution, the former 
being quite conspicuous. The dotted lines are positioned at 
log Knax = 2.55 and log V;„ax = 1.675 to guide the eye. Fig. |] 



log Dlinear 

Fig. 5. The inverse Tully-Fisher diagram for the sample used 
in the analysis. The solid line refers to a linear regression of 
a' — 0.576 and b' = 1.256. The dashed lines give the forced 
solutions with a' = 0.50 for Hubble types 1 with b' = 1.448 
and 8 with b' = 1.209 The dotted fines at log Kiax = 2.55 and 
log Vniax = 1.675 are intended to guide the eye. At least the 
upper cut-off is quite conspicuous. 
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Fig. 6. The differential behaviour of (logiJo) as a function of 
the normalized distances. The inverse parameters were a' = 0.5 
and b' = 1.450. 



hints that the slope a' — 0.5 adopted in Sect. 3 may not be im- 
peccable and thus questions the validity of the "naive" deriva- 
tion of Hq at least in the case of the diameter sample. 

In the case of diameter samples, Teerikorpi et al. ( 1999| ) 



discussed how the cut-offs should demonstrate themselves in a 
logiJo vs. logdnorm diagram, where log(i„orm = log£'25 + 
logdkin, which in fact is the log of Dunear next to a constant. 
We call dnorm "normalized" in analogy to the method of nor- 
malized distances, where the kinematical distances were nor- 
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log ^Ceph 

Fig. 7. A straightforward linear regression applied to the cali- 
brator sample yielding a' = 0.749 and b' = 1.101. 



Fig. 9. The differential (log_ffo) vs. /Xnorm 
no indication of a selection in logKnax- 
(a' = —0.10) appears to be incorrect. 



A 
o 
X 

cn 
_o 

V 



a' = 0.749, b' = 



1099 



0) in 



1 .5 



2.5 



loq d 



Fig. 8. As Fig. but now the parameters a' = 0.749 and 
b' = 1.101 were used. One can see how the downward ten- 
dency between log dnorm ^ 1-45 and log dnorm ~ 2 has disap- 



peared. Also cf. Fig. 2 in Teerikorpi et al. ( |1999[ ) 



malized in order to reveal the underlying bias. That is exactly 
what is done also here. 

Consider the differential behaviour of (log Hq) as a func- 
tion of the normalized distance. Differential average (logiJo) 
was calculated as follows. The abscissa was divided into inter- 
vals of 0.01 starting at minimum logdnorm in the sample. If a 
bin contained at least 5 galaxies the average was calculated. In 
Fig. ^. the inverse parameters a' = 0.5 and 6' = 1.450 were 
used. It is seen that around log dn ^ 2 the values of log Hq 
have a turning point as well as at logdn ~ 1.45. The most 
striking feature is - however - the general decreasing tendency 
of logiJo between these two points. Now, according to ET97, 
a downward tendency of log Hq as a function of distance cor- 



diagram. One finds 
The adopted slope 



responds to A/ A' > 1, i.e. the adopted slope A is too shallow 
(A' is the relevant slope). 

Closer inspection of Fig. |l] shows that a steeper slope might 
provide a better fit to the calibrator sample. One is thus tempted 
to ask what happens if one adopts for the field sample the slope 
giving the best fit to the calibrator sample. As such solution 
we adopt the straightforward linear regression yielding a' = 
0.749 and b' = 1.101 shown in Fig. ^ It is interesting to note 
that when these parameters are used the downward tendency 
between logdnorm ^ 1-45 and logdnorm ^ 2 disappears as 
can be seen in Fig. || From hereon we refer to this interval 
as the "unbiased inverse plateau". The value of logiJo in this 
plateau is still rather high. 

In the case of the magnitude sample we study the behaviour 
of the differential average (logi/o) as a function of a "normal- 
ized" distance modulus: 



Mil 



- 51ogdkin. 



(22) 



The fJ-noim axis was divided into intervals of 0.05 and again, if 
in a bin is more than five points the average is calculated. As 
suspected in the view of Fig. 0., Fig. H reveals no significant 
indications of a selection in log Kiiax- The points follow quite 
well the straight line also shown. The line however is tilted 
telling us that the input slope a' — —0.10 may not be the rele- 
vant one. 

As already noted the type corrections may have some influ- 
ence on the slopes. In the next section we derive the appropriate 
type corrections for the zero-points using galaxies residing in 
the unbiased plateau (log dnorm G [1.45, 2.0]) for the diameter 
sample and for the whole magnitude sample and rederive the 
slopes. 
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Ab'{T) 


a' = 0.54 


a' = -0.115 


Ab\l) 


0.125 


0.110 


A6'(2) 


0.156 


0.124 


A6'(3) 


0.129 


0.096 


A6'(4) 


0.095 


0.058 


A6'(5) 


0.069 


0.030 


Afe'(6) 


0.0 


0.0 


Ab'{7) 


-0.054 


-0.042 


A6'(8) 


-0.118 


-0.075 



Table 1. The type corrections required for the relevant slopes 
a' — 0.54 for the unbiased diameter sample and a' = —0.115 
for the magnitude sample. 
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I ' 
a' = 0.54 
b' = 1.3252 



5. Type corrections and the value of Ho 

The zero-points needed for the type corrections are calculated 
using Eq. 2l]. It was pointed out in Sect. 2.2 that log Ho should 
run horizontally in order to find an unbiased estimate for Hq. 
In this section we look for such an slope. Because the type- 
corrections depend on the adopted slope, this fine-tuning of the 
slope must be carried out in an iterative manner This process 
consists of finding the type corrections A&'(T) for each test 
slope a'. Corrections are made for both the field and calibrator 
samples. The process is repeated until a horizontal (log Hq) run 
is found. 

Consider first the diameter sample. When the criteria for 
the unbiased inverse plateau were imposed on the sample, 2142 
galaxies were left. For this subsample the iteration yielded 
a' = 0.54 (the straight line in Fig. |l^ is the least squares fit 
with a slope 0.003) and when the corresponding type correc- 
tions given in Table 1 were applied to the calibrator sample and 
the slope forced to it one found b'^^i{6) = 1.325. The result 
is shown in Fig. |l^. The given inverse parameters predict an 
average (logi/o) = 1-897 (or Hq = 78.9 kms^i Mpc~i). 

We treated the magnitude sample of 1713 galaxies in a sim- 
ilar fashion. The resulting best fit is shown in Fig. 11 The 
relevant slope is a' — —0.115 (the least squares fit yields a 
slope 0.0004). The corresponding type corrections are given 
in Table 1. The forced calibration gives b'^^^{6) — —0.235. 
From this sample we find an average {log Ho) — 1.869 (or 
Hq = 72.4 kms~^ Mpc~^). In both cases the inverse esti- 
mates for the Hubble constant (Ho « 80 for the diameter rela- 
tion and Hq « 70 for the magnitude relation) are considerably 
larger than the corresponding estimates using the direct TuUy- 
Fisher relation (Hq w 55). 

6. Ho corrected for a calibrator selection bias 

The values of Hq from the direct and inverse relations still 
disagree even after we have taken into account the selection 
in log Knax, made the type corrections and used the relevant 
slope. There is - however - a serious possibility left to ex- 
plain the discrepancy. The calibrator sample used may not meet 
the theoretical requirements of the inverse relation. In order 
to transform the relative distance scale into an absolute one a 



Fig. 10. The differential (log Hq) as a function of the log of 
normalized distance log dnoim for the plateau galaxies with the 
adopted relation log Vmax — 0.54 log D + 1.325. The solid line 
is the average logi?o ~ 1.897. 
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Fig. 11. The differential (logiJo) 
normalized distance modulus /i„, 
with the adopted relation log Vma; 
solid line is the average log Hq = 



as a function of the log of 
,rin for the plateau galaxies 
; = -0.115M- 0.235. The 
1.869. 



properly chosen sample of calibrating galaxies is needed. What 
does "properly chosen" mean? Consider first the direct relation 
for which it is essential to possess a calibrator sample, which 
is volume-limited for each pcai- This means that for a pcai one 
has Xcax which is drawn from the complete part of the gaus- 
sian distribution function G{X\Xp,(Jx ), where the average 
Xp = ap + b. If ax is constant for all p and the direct slope 
a has been correctly derived from the unbiased field sample, it 
will, when forced onto the calibrator sample, bring about the 
correct caUbrating zero-point. 

As regards the calibration of the inverse relation the sample 
mentioned above does not necessarily guarantee a successful 
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calibration. As pointed out by Teerikorpi et al. (1999) though 
the calibrator sample is complete in the direct sense nothing 
has been said about how the Pcai's relate to the correspond- 
ing cosmic distribution of p's from which the field sample was 
drawn, (p)cai should reflect the cosmic average pq. If not, the 
relevant field slope when forced to the calibrator sample will 



bring about a biased estimate for Hq . Teerikorpi ( 1 990 ) akeady 



recognized that this could be a serious problem. He studied, 
however, clusters of galaxies where a nearby (calibrator) clus- 
ter obe ys a different slope than a distant cluster. Teerikorpi et al. 
( 1999 ) developed the ideas further and showed how this prob- 
lem may be met also when using field galaxies. The mentioned 
bias when using the relevant slope can be corrected for but is a 
rather complicated task. For the theoretical background of the 
"calibrator selection bias" consult Teerikorpi et al. (1999). 

One may - as pointed out by Teerikorpi et al. (1999) - use 
instead of the relevant slope the calibrator slope which also pre- 
dicts a biased estimate for Hq but which can be corrected for in 
a rather straightforward manner. For the diameter relation the 
average correction term reads as 



AlogiJo = (3-a) InlOcr 



1 



(23) 



where ajj is the dispersion of the log linear diameter 
log Z^iincar and a gives the radial number density gradient : 
a = corresponds to a strictly homogeneous distribution of 
galaxies. For magnitudes the correction term follows from (cf. 
Teerikorpi |1990|) 



AlogiJo = 0.2 



"-cal 

a' 



1 



X i(M) - Ma) 



(24) 



Because — A/q simply reflects the classical Malmquist bias 
one finds: 



Ai rr (3 -a) In 10 2 
Alogi/o = ^ r '^M X 0.2 



-'cal 



1 



(25) 



Note that one may use the calibrator slope and consequently 
the correction formulas irrespective of the nature of the cali- 
brator sample (Teerikorpi et al. 1999 1. If the calibrator sam- 



ple would meet the requirement mentioned, the value corrected 
or H should equal values obtained from the rele- 



with Eqs. 23 



vant slopes. Furthermore, our analysis carried out so far would 
have yielded an unbiased estimate for Hq and thus the prob- 
lems would be in the direct analysis. However, if the require- 
ment is not met one should prefer the corrective method using 
the calibrator slope. 

6.1. Is the calibrator sample representative? 

Is the calibrator bias present in our case? Recall that the cali- 
brators used were sampled from the nearby field to have high 
quality distance moduli mostly from the HST Cepheid mea- 
surements. This means that we have no a priori guarantee that 
the calibrator sample used will meet the criterium required. We 
compare the type-corrected diameter and magnitude samples 
with the calibrator sample. Note that for the diameter sample 
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D 
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log V^ax 

Fig. 12. Histogram of the log T^nax values and the individual 
calibrators (labelled with stars). The vertical solid line gives the 
median of the plateau Med(log V^i^ = 2-10 and the dot- 
ted line gives the median of the calibrators Med(log V^-^^^) — 
2.11. 
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Fig. 13. The Kolmogorov-Smirnov test for the diameter sam- 
ple. Pay attention to the rather remarkable similarity between 
the cumulative distribution functions (cdfs). 



we use only galaxies residing in the unbiased inverse plateau 
(i.e. the small selection effect in log K„ax has been eliminated). 
In Fig. |l^ we show the histogram of the log Vmax values for 
the diameter sample and the individual calibrators (labelled as 
stars). The vertical solid line gives the median of the plateau 
Mcd(log V„p!^'^'°''") = 2.10 and the dotted line gives the me- 
dian of the calibrators Med (log V,,';!'^'') = 2.11. In the case of 
magnitudes both the field and calibrator sample have the same 
median (2.14). The average values for the diameter case were 
(logViP|f^'°^") = 2.09 and {\ogV,^^^^) = 2.06, and for the 
magnitude case {\ogV^^^) = 2.12 and {logV;;^^^) = 2.08. 
Both the diameter and the magnitude field samples were sub- 
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Fig. 14. The Kolmogorov-Smirnov test for the magnitude sam- 
ple. Again the cdfs are quite similar. 



jected to strict limits, which means that both inevitably suffer 
from the classical Malmquist bias. In order to have a repre- 
sentative calibrator sample in the sense described, we would 
have expected a clear difference between the field and calibra- 
tor samples. That the statistics are very close to each other lends 
credence to the assumption that the calibrator selection bias is 
present. 

We also made tests using the Kolmogorov-Smirnov statis- 
tics (Figs. [l3|and[l4|). In this test a low significance level should 
be considered as counterevidence for a hypothesis that two 
samples rise from the same underlying distribution. We found 
relatively high significance levels (0.89 for the diameter sam- 
ple and 0.3 for the magnitude sample). Neither these findings 
corroborate the hypothesis that the calibrator sample is drawn 
from the cosmic distribution and hence the use of Eqs. ^ or ^ 
is warranted. 

6.2. The dispersion in log Diincar 

In order to find a working value for the dispersion in 
logDiinoai , we first consider the classical Spaenhauer diagram 
(cf. Sandage 1994a, 1994b| ). In the Spaenhauer diagram one 
studies the behaviour of X as a function of the redshift. If the 
observed redshift could be translated into the corresponding 
cosmological distance, then X inferred from x and the redshift 
would genuinely reflect the true size of a galaxy. 

In practice, the observed redshift cannot be considered as 
a direct indicator of the cosmological distance because of the 
inhomogeneity of the Local Universe. Peculiar motions should 
also be considered. Thus the inferred X suffers from uncer- 
tainties in the underlying kinematical model. The Spaenhauer 
diagram as a diagnostics for the distribution function is always 
constrained by our knowledge of the form of the true velocity- 
distance law. 

Because the normalized distance (cf. Sect. 3.) is propor- 
tional to the linear diameter we construct the Spaenhauer dia- 




Fig. 15. A classical Spaenhauer diagram for normalized dis- 
tances vs. kinematical distances with a presumed dispersion 
ax = 0.28. 
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Fig. 16. Comparison between average values of (log dnorm) for 
different kinematical distances and the theoretical prediction 
calculated from Eq. |6|with X'^ = 1.37 and ax = 0.28. 

gram as log dnorm vs. log dkin thus avoiding the uncertainties in 
the absolute distance scale. The problems with relative distance 
scale are - of course - still present. The fit shown in Fig. [IJ 
is not unacceptable. The dispersion used was ax — 0.28, a 
value inferred from the disper sion in absolute B-band magni- 
tudes aM = 1.4 (Fouque et al. 1990 ) based on the expectation 
that the dispersion in log linear diameter should be one fifth of 
that of absolute magnitudes. 

We also looked how the average values (logdnorm) at dif- 
ferent kinematical distances compare to the theoretical predic- 
tion which, in a strictly limited sample of X's, at each log dis- 
tance is formally expressed as 



= X* + 



2ax cxp[-(Xiin,-Xo*)V(2a|,)] 
V2ti erfc[(Xiin, - X*)/{V2ax)\ 



(26) 
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Fig. 17. A least squares fit the type corrected calibrator sample 
yielding a' = 0.73. The type correction was based on a' = 
0.54. 
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Fig. 18. A least squares fit the type corrected calibrator sample 
yielding a' = —0.147. The type correction was based on a' = 
-0.115. 



Here X refers to logd„orm- The curve in Fig. 16 is based on 
Xq = 1.37 and ax = 0.28. The averages from the data are 
shown as bullets. The data points follow the theoretical predic- 
tion reasonably well. 

6.3. Corrections and the value of Hq 

Consider a strictly homogeneous universe, i.e. a = 0. In Eqs. 
^ and ^ one needs values for slope a^. Least squares fit to 
the type-corrected calibrator sample yields a'^ = 0.73 for the 
diameter relation and a'^ — —0.147 for the magnitude relation, 
(cf. Figs. [T^ and |l8|). These slopes correspond to diameter zero- 
point 5^(6) — 1.066 ± 0.103 and to magnitude zero-point h'^ = 
—0.879 ± 0.131 The biased estimates for average logiJo are 
(logi/o) = 1.910 ± 0.188 for the diameters and (logi/o) = 



1.876 ± 0.176 for the magnitudes. For the zero-points and the 
averages we have given the la standard deviations. The mean 
error in the averages is estimated from 



e(iogHo> 



'log Ho 



Nr. 



gal 



(27) 



27 



where as' = at, / a^^j 
for magnitudes. The use of Eq 
dispersion in b' does not correlate with the dispersion logifo- 
With the given slopes and dispersions we find: 



for diameters and as' = 0.2ab'/a'^^^ 
is acceptable because the 



- (log Ho) 

- (logi/o) 



1.910 ± 0.037 for the diameters 
1.876 ± 0.046 for the magnitudes. 




Eq. ^ predicts an average correction term for the slopes 
< = 0.73 and a' = 0.54 together with ax = 0.28 A log Hq = 
0.191. and Eq. || with 4 = -0.147,a' = -0.115 and an = 
1.4 AlogiJo — 0.151. When applied to the above values we 
get the corrected, unbiased estimates 

- (log Ho) = 1.719 ± 0.037 for the diameters 

- (log Ho) = 1.725 ± 0.046 for the magnitudes. 

These values translate into Hubble constants 

- Ho = 52^4 kms~^ Mpc"^ for the inverse diameter B- 
band Tully-Fisher relation, and 

- Ho = 53^5 kms^^ Mpc^^ for the inverse magnitude B- 
band Tully-Fisher relation. 

These corrected values are in good concordance with each 
other as well as with the estimates established from the direct 
diameter Tully-Fisher relation (Theureau et al. 1997b). Note 



that the errors in the magnitude relation are slightly larger than 
in the diameter relation. This is expected because for the diam- 
eter relation we possess more galaxies. The error is however 
mainly governed by the uncertainty in the calibrated zero-point. 
This is expected because though the dispersion in inverse rela- 
tion as such is large it is compensated by the number galaxies 
available. 

Finally, how significant an error do the correction formu- 
lae induce? We suspect the error to mainly depend on a. The 
correction above was based on the assu mption of homogeneity 
(i.e. a = 0). Recently Teerikorpi et al. ( 1998 ) found evidence 
that the average density radially decreases around us {a w 0.8) 
confir ming the more general (fractal) analysis by Di Nella et al. 
( |1994 ). Using this value of a we find A log = 0. 140 for the 
diameters and A log Ho = 0.111 for the magnitudes yielding 

- (log Ho) = 1.770 ± 0.037 for the diameters 

- (log Ho) = 1.765 ± 0.046 for the magnitudes. 

In terms of the Hubble constant we find 



- Ho = 59_4kms Mpc for the inverse diameter B- 
band Tully-Fisher relation, and 

- Ho ^ b^tl kms-i Mpc"! for the inverse magnitude B- 
band Tully-Fisher relation. 
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7. Summary 

In the present paper we have examined how to apply the inverse 
Tully-Fisher relation to the problem of determining the value 
of the Hubble constant, Hq, in the practical context of the large 
galaxy sample KLUN. We found out that the implementation 
of the inverse relation is not as simple task as one might expect 
from the general considerations (in particular the quite famous 
result of the unbiased nature of the relation). We summarize 
our main results as follows. 



1. A straightforward application of the inverse relation con- 
sists of finding the average Hubble ratio for each kinemat- 
ical distance and tranforming the relative distance into an 
absolute one through calibration. The 15 calibrator galaxies 
used were drawn from the field with cepheid distance mod- 
uli obtained mostly from the HST observations. The inverse 
diameter relation predicted Hf) w 80 kms~^ Mpc~^ and 
the magnitude relation predicted Hq w 70 kms~^ Mpc~^ 
The diameter value for Hq is about 50 percent and the mag- 
nitude value about 30 percent larger than those obtained 
from the direct relation (cf. Theureau et al. 1997b). 



3. 



We examined whether this discrepancy could be resolved 
in terms of some selection effect in log Vmax and the type 
dependence of the zero-points on the Hubble type. One ex- 
pects these to have some influence on the derived value of 
Hq. Only a minuscule effect was observed. 
There is - however - a new kind of bias involved: if the 
log Vniax-distribution of the calibrators does not reflect the 
cosmic distribution of the field sample and the relevant 
slope for the field galaxies differs from the calibrator slope 
the average value of log Hq will be biased if the relevant 



slope is used (Teerikorpi et al. 1999) 



4. We showed for the unbiased inverse plateau galaxies i.e. a 
sample without galaxies probably suffering from selection 
in log Vmax, that the calibrators and the field sample obey 
different inverse diameter slopes, namely a^^j — 0.73 and 
a' = 0.54, Also, the magnitude slopes differed from each 
other (a[,^j = —0.147 and a' = —0.115). For the diameter 
relation we were able to use 2142 galaxies and for the mag- 
nitude relation 1713 galaxies. These sizes are significant. 

5. We also found evidence that the calibrator sample does 
not follow the cosmic distribution of log Vmax for the field 
galaxies. This means that if the relevant slopes are used a 
too large value for Hq is found. Formally, this calibrator 
selection bias could be corrected for but is a complicated 
task. 

6. One may use instead of the relevant slope the calibrator 
slope which also brings about a biased value of Hq. Now, 
however, the correction for the bias is an easy task. Further- 
more, this approach can be used irrespective of the nature 
of the calibrator sample and should yield an unbiased esti- 
mate for Hq. 

7. When we adopted this line of approach we found 

- -ffo = 52^4 km Mpc~^ for the inverse diameter B- 
band Tully-Fisher relation, and 



- Hq = 53^5 kms~^ Mpc~^ for the inverse magnitude 
B-band Tully-Fisher relation 

for a strictly homogeneous distribution of galaxies (a = 0) 
and 

- Hq — 59^4 km s^^ Mpc~^ for the inverse diameter B- 
band Tully-Fisher relation, and 

- Hq = 58^5 kms~^ Mpc~^ for the inverse magnitude 
B-band Tully-Fisher relation. 

for a decreasing radial density gradient (a = 0.8). 



These values are in good concordance with each other as 
well as with the values established from the correspon ding di - 
rect Tully-Fisher relations derived by Theureau et al. (1997b), 
who gave a strong case for the long cosmological distance scale 
consistently supported by Sandage and his collaborators. Our 
analysis also establishes a case supporting such a scale. It is 
worth noting that this is the first time when the inverse Tully- 
Fisher relation clearly lends credence to small values of the 
Hubble constant Hq. 
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Appendix A: The relevant slope and an unbiased Hq 

In this appendix we in simple manner demonstrate how the rel- 
evant slope introduced in Sect. 2.2 indeed is the slope to be 
used. Consider 



logi/o log Vcor - l0gi?iTF, 



(A.l) 



where the velocity corrected for the peculiar motions, Vcor, de- 
pends on the relative kinematical distance scale as 

log Vcor = log Ci + log dkin (A.2) 

and the inverse Tully-Fisher distance in Mpc is^ 

log i?iTF ^Ap + Seal -x + f3 (A.3) 

The constant Ci was defined by Eq. [ll| and can be decomposed 
into logCi ~ logi?o + logC2. Hq is the true value of the 
Hubble constant and C2 transforms the relative distance scale 
into the absolute one: log_Rki n = logdkin + logC2. Because 
-^kin — log i?kin + X — (3 Eq. Al rcads: 



log Hq - log H* = Xkin -Ap-B, 



cal • 



(A.4) 



Consider now a subsample of galaxies at a constant Xq. By 
realizing that Xkin = Xq + [B' — i?in), where B' gives the 
true distance scale and i?in depends on the adopted distance 



^ The numerical constant /3 = 1.536274 connects a; in 0.1 arcsecs, 
X in kpc and i?iTF in Mpc 
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scale (based on the input Hq), and by taking the average over 
Xo Eq. |A.4| yields 



(log Ho)x^ - log H* = Xkin - A{p)x^ - Seal. (A.5) 

The use of B' is based on two presumptions, namely that the 
underlying kinematical model indeed brings about the correct 
relative distance scale and that the adopted value for Ci gen- 
uinely reflects the true absolute distance scale. If the adopted 
slope is the relevant one we find using Eq. |l9| A{p)xa — 
Xo — -Bin and 

{log Ho) - logH* = {B' - Bin) - (Seal - ^in). (A.6) 

As a final result we find 



{\ogHo)x^-\ogH* = {b'^ 



cal 



6;,,J/<. (A.7) 



Because Eq. |A.7| is valid for each Xo, the use of the relevant 
slope necessarily guarantees a horizontal run for (log iJo) as a 
function of Xq. 

Appendix B: Note on a theoretical diameter slope 
a' ~ 0.75 



Theureau et al. ( |1997a| ) presented theoretical arguments which 
supported the inverse slope a' — 0.5 being derived from the 
field galaxies. Consider a pure rotating disk (the Hubble type 
8). The square of the rotational velocity measured at the radius 
'"max at which the rotation has its maximum is directly propor- 
tional to the mass within r-^^^, which in turn is proportional to 
the square of r^ax- Hence, log Kiax « 0.5 logrmax- By adding 
a bulge with a mass-to-luminosity ratio differing from that of 
the disk, and a dark halo with mass proportional to the lumi- 
nous mass, one can as a first approximation understand the de- 
pendence of the zero-point of the inverse relation on the Hubble 
type. 

However, the present study seems to require that the the- 
oretical slope a' is closer to 0.75 rather than 0.5. The ques- 
tion arises whether the simple model used by Theureau et al. 



( 1997a) could in some natural way be revised in order to pro- 



duce a steeper slope. In fact, the model assumed that for each 
Hubble type the mass-to-luminosity ratio M/L is constant in 
galaxies of different sizes (luminosities). If one allows M/L 
to depend on luminosity, the slope a' will differ from 0.5. Es- 
pecially, if M/L oc L^-^^, one may show that the model pre- 
dicts the inverse slope a' = 0.75. The required luminosity de- 
pendence of M/ L is interestingly similar to that of the funda- 
mental plane for elliptical galaxies and bulges (Burstein et al. 



1997[ ). The questions of the slope, the mass-to-luminosity ratio 



and type-dependence will be investigated elsewhere by Hanski 
& Teerikorpi (1999, in preparation). 



Appendix C: How to explain a'^^^ < a"! 

Among other things ET97 discussed how a gaussian measure- 
ment error (j^ in apparent diameters yields a too large value 
for Hq. How does the combination of cut-offs in logV^^ax - 
distribution and this bias affect the slope? We examined this 
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Fig. C.l. A synthetic Virgo supercluster subjected to a bias 
caused by measurement errors in apparent diameters and upper 
and lower cut-offs in log Vinax- 

problem by using a synthetic Virgo Supercluster (cf. Ekholm 
|1996 ). As a luminosity function we chose 



logD = 0.28 X G{0,1) 
and as the inverse relation 



1.2 



(C.l) 



logKnax = 0.11 X G(0,1) +a;iogi:' + 0.9, (C.2) 

where G(0, 1) refers to a normalized gaussian random variable. 
As the "true" inverse slope we used a(. — 0.75. The other nu- 
merical values were adjusted in order to have a superficial re- 
semblance with Fig. ^. We first subjected the synthetic sample 
to the upper and lower cut-offs in log Vniax given in Sect. 4. 
The resulting slope was a' = 0.692. A dispersion of ax = 0.05 
yielded a' = 0.642 and ax = 0.1 a' = 0.559. The inverse 
Tully-Fisher diagram for the latter case is shown in Fig. |C.lt 
Though the model for the errors is rather simplistic this exper- 
iment shows a natural way of flattening the observed slope a' 
with respect to the input slope a[ . 
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